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CONVECTION OF A PATTERN OF VORTICITY THROUGH A SHOCK WAVE ' 

By H. S. Kibneb 


SUMMARY 

An arbitrary weak spatial distribution of vorticity can be 
represented in terms of plane sinusoidal shear waves of all 
orientations and wave lengths {Fourier integral). The analysis 
treats the passage of a single representative weak shear wave 
through a plane shock and shows refraction and modification of 
the shear wave with simultaneous generation of an acoustically 
intense sound wave. Applications to turbulence arid to noise 
in supersonic wind tunnels are indicated. 

INTRODUCTION 

Turbulence such as the residual small oddjung motion in a 
wind-tunnel stream will gradually decay as it is carried along. 
The decay process has been the subject of much study in the 
face of formidable difficulties. The random character of the 
motions has been successfully handled by the methods of 
statistics; even with these methods, however, the non- 
linearity of the equations governing the intermi.xing processes 
has severely limited the progress attainable without sim- 
plifying assumptions. 

On the other hand, for relatively sudden changes in 
turbulence, such as occur when it passes through a wire- 
mesh damping screen, the decay may be negligible and the 
changes may follow linear laws. The linearity is assured if 
the turbulence constitutes a sufficiently small perturbation 
of the main stream. Recently it has been found that the 
problem of such linear changes could be solved completely 
by a specialized adaptation of the spectrum concept of the 
statistical theory of turbulence. 

Several of these linear processes have been treated in this 
manner; the damping-screen problem (ref. 1), the passage of 
turbulence through a sudden wind-tunnel contraction (ref. 2), 
and the passage of turbulence tlu’ough a series of screens 
followed by a sudden contraction (ref. 3). A basic tech- 
nique for such problems has been evolved in these papers. 

The present paper is motivated by another problem of 
the same linear character, namely, the convection of weak 
turbulence through a shock wave. Among other circum- 
stances, this problem arises in the interpretation of measure- 
ments with a hot-wire anemometer in a supersonic stream, 
because a detached bow wave stands ahead of the wire.* 

1 Supersedes NACA TN 2864, "Convection of a Pattern of Vorticity Through a Shock 
Wave” by H. S. Rlbner, 1953. 

2 A simple interpretation for all but very small eddies comparable with the scale of the bow 
wave is, however, available In the work of Kovdsznay (ref. 4). 


Such a curved shock is not attractive for theoretical analysis, 
hut it is not difficult to replace it with an extended plane 
shock by use of auxiliary means; attention can thus be 
limited to the convection of turbulence through a plane 
shock. 

The coneeptual basis for the treatment of these linear 
problems is as follows: An arbitrary weak spatial distribu- 
tion of vorticity — and hence a weak turbulent velocity 
field — can be represented as a superposition or spectrum of 
plane sinusoidal shear waves distributed among all orienta- 
tions and wave lengths. This is a physical interpretation 
of the mathematical formulation as a Fourier integral; ^ the 
individual shear waves may be identified as Fourier or 
spectrum components. When the turbulence wave pattern 
is convected through a screen or through a shock wave, the 
individual waves are altered without mutual interference if 
the waves are suitably weak. Thus the modified field 
downstream of the screen or shock can be obtained, in 
principle, by superposition of the modified individual waves. 
In practice the description of the detailed spatial distribu- 
tion of velocity, either initially or finally, is hopeless; the 
initial wave distribution is known onl_y statistically (e. g., 
the phase angles are unknown), and statistical changes only 
can be calculated. In either case the analysis of the be- 
havior of a representative single wave constitutes a pre- 
requisite to the determination of the changes in the weak 
turbulent field. 

In the present paper such an analysis is carried out for a 
single shear wave, of arbitrary inclination, convected through 
a plane shock. There remains the task of calculating there- 
from the changes in the statistical properties of a weak turbu- 
lent field convected through a shock. Suitable procedures 
have been developed in references 1 and 2 ; for their applica- 
tion the present results, which are formulated in two dimen- 
sions, must be reexpressed to bring out the spatial inclination 
of the wave in three dimensions. The procedures will re- 
quire some modification before the noise field generated by 
the interaction can be treated. 

3 The velocity field so represented may be either rotational or irrotational within the speci- 
fied region, even though the “building blocks,” the shear waves, are rotational. In case an 
irrotational field is represented, the vorticity of these shear waves, but not the velocity, 
mutually cancels within the specified region (which may be multiply connected), leaving a 
distribution of vorticity in the external space. The irrotational flow may be regarded as 
Induced by this external vorticity. 

These remarks all refer to a velocity field satisfying the incompressible continuity equation: 
a smalI-i>erturbation field of vorticity in fluid at rest, or convected by a main stream, will 
fulfill this condition. 
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Ficure 1 . — Convection of plane oblique sinusoidal shear wave through 
shock: original unsteady-flow problem. 

This single-wave problem is also treated in reference 5. 
The analyses bear little resemblance; In that work a refer- 
ence frame is used in which the flow is unsteady, whereas 
herein a frame is used in which the flow is steady. Sound 
waves are likewise treated in the work cited. 

The outline of the present analysis is as follows: The prob- 
lem is posed as the calculation of the flow field behind a plane 
normal shock wave due to the convection through the shock 
of an inclined plane sinusoidal shear wave; the shear wave is 
specified to be weak to ensure small perturbations to the 
mean flow. This problem, for which the flow is unsteady in 
time, is converted into an equivalent steady-flow problem 
by transformation to a moving frame of reference. In this 
frame the normal shock is replaced by an equivalent oblique 
shock. 

The analysis is now formulated as a boundary-value prob- 
lem for the flow in the region downstream of the shock: The 
governing partial differential equation for this small-pertur- 
bation rotational flow is derived (extension of Seam’ work, 
ref. 6) ; boundary conditions on the velocity components just 
behind the shock are obtained from the oblique-shock rela- 
tions; and finally the rotation term in the governing equation 
is evaluated in terms of gradients of entropy and total 
enthalpy, with use of the entropy changes across the shock. 
The initially unknown perturbation of the form of the shock 
wave is taken into account in the boundary conditions and 
rotation term by assuming it to be sinusoidal with initially 
undetermined amplitude and phase. 

The velocit}' W (all symbols are defined in appendix A) 
downstream of this equivalent oblique shock may be either 
subsonic or supersonic depending on the inclination of the 
initial plane shear wave; separate solutions of the boundary- 
value problem are worked out for the two markedly different 
cases. The horizontal shear wave — which is a simple 
special case for subsonic W — is given a separate treatment. 

The analysis of the velocity field downstream of the shock 
is followed by an account of the associated pressure, density, 
and temperature fields there and of the distortion of the ini- 
tially plane shock. Finallj", the acoustic level of the (fluc- 
tuating) pressure field is worked out in approximate fashion 
for an example applied to a supersonic wind tunnel: A par- 
ticular initial intensity of turbulence is assumed and con- 
sidered as being concentrated in a single shear wave rather 



FiotiRB 2. — Transformation to equivalent steady-flow problem by 
superposition of velocity V. 


than distributed throughout a continuous spectrum. The 
calculation amounts to an estimate of the noise level generated 
by the passage of a specified level of turbulence through a 
shock wave. 

This investigation was conducted at the NACA Lewis 
laboratory. 

FORMULATION OF BOUNDARY-VALUE PROBLEM 

The unsteady-flow problem. — The inclined plane sinu- 
soidal shear wave is shown schematically in figure 1. The 
flow is viewed in a plane perpendicular to the shock and to 
the wave fronts. The wave is supposed to be convected 
downstream by the main stream with velocity so that it 
passes through the normal shock. The passage through the 
shock is evidently an unsteady process, since the intercepts 
of the inclined lines (the nodes of the sine wave) move down- 
ward along the shock front; it will be shown that a disturbance 
ripple moves along the shock with the same speed V. 

In the general case of a plane oblique sinusoidal shear 
wave there will also exist a perturbation velocity component 
normal to the plane of the figure. Now the ripples in the 
shock front will be two dimensional, and the shock with the 
ripples will still be everywhere perpendicular to the plane of 
the figure. Thus, the normal velocity component will be 
parallel everywhere to the shock and will be unaffected as the 
shear wave passes through; the component will have no other 
effect. Its invariance established, this normal velocity com- 
ponent will be omitted from the analysis. 

The equivalent steady-flow problem. — If an observer 
moves downward along the shock with a speed V, relative 
to him the flow will have an apparent upward velocity com- 
ponent V. This scheme of things is shown in figure 2. In 
particular, V has been chosen so that the resultant stream 
velocity (relative to the moving observer) is alined with the 
velocity in the disturbance wave; that is, V— Ua tan d. The 
observer then sees what appears to be a steady sinusoidal 
shear flow passing through an oblique shock wave. Thus, by 
the proper choice of a system of moving axes the original 
unsteady-flow problem has been converted into an equivalent 
steady-flow problem. 
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FinuBE 3. — Symbols and coordinate axes. 


Governing partial differential equation for rotational 
flow. — The task of the analysis is to calculate the flow field 
on the downstream side produced by the passage of the sinus- 
oidal shear flow through the equivalent oblique shock. It is 
to be expected that the shock will be perturbed from its mean 
plane and will, in fact, develop a corrugated appearance. 
Because of these corrugations, vorticity (rotation) will be 
introduced into the downstream flow. This vorticity and all 
the downstream velocity perturbations will be weak com- 
pared with the stream velocity because the original disturb- 
ance wave has been assumed weak. Thus, a small-perturba- 
tion, or linearized, treatment of the flow field is permissible. 

In reference 6 the governing partial differential equation 
for small-perturbation compressible rotational flow has been 
derived for isoenergetic flow, that is, for flow of constant 
stagnation enthalpy. However, the shear wave under con- 
sideration possesses variable stagnation enthalpy; that is, 
pressure, density, and temperature are constant upstream of 
the shock, but the velocity varies. It has been necessary, 
therefore, to obtain a more general governing equation that 
applies when both entropy and stagnation enthalpy are 
variable. The derivation is given in appendix B. 

This governing equation is expressed in terms of coordi- 
nates ^ and 7}, I being the distance in the main stream direction 
and -q the distance perpendicular thereto. The equation 
reads 

( 1 ) 

where W is the stream velocity in the transformed problem, 
W is the corresponding Mach number, H is the stagnation 
enthalpy, s is the entropy, T is the temperature, is the 
vorticity, and ^ is a perturbation stream function such that 
^7j = w=perturbation velocity in J direction 

— — = = perturbation velocity in direction 



Unperturbed 

Perturbed 

y 



Fioure 4. — Geometrical relations across shock, with and without 
perturbation <r in shock angle. 


axes. In this figure 14^ is the resultant stream velocity down- 
stream of the shock (in the moving frame of reference), and 
the f and rj axes are indicated. The final flow pattern depends 
crucially on whether W is subsonic or supersonic ; the criterion 
depends, in turn, on the Mach number corresponding to Ua 
and on the wave inclination 6. 

Boundary conditions. — The boundary conditions just 
downstream of the shock will now be obtained by application 
of the shock-wave relations. 

By geometiy (fig. 4) the stream velocity components 
normal and tangential to the undisturbed shock arc, respec- 
tively, 

Ua=Wa cos 6 


V= Wa sin 6 


The shear wave will provide directly a perturbation Wa to Wa 
and will cause indirectly a perturbation (r{y) to the shock- 
wave angle, of initially undetermined magnitude. The 
effect of er is equivalent to an increment in 6. The asso- 
ciated perturbations to Ua and V are found by obtaining 
their respective differentials and replacing (IWa by Wa and dd 
by ff therein; the results are 


dUA=WA cos d — (tWa sin 8 
dV=WA sin 8-\-(tWa cos 8. 


( 3 ) 


The corresponding change in normal velocity U downstream 
of the shock is obtained from the normal-shock relation 


(The stream function is defined differently in ref. 6, as it 
involves an entropy term.) 

For application of equation (1) in the present problem 
reference shoidd he made to figure 3 for the dii’ection of the 


Ua 

U 
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By logarithmic differentiation and use of the fact that the 

upstream temperature is constant (^whence 

\ Ua ^ ^ ' 

there is finally obtained 


dU dUA(. 
U~ UaV 


-2 


7—1 

7+1 



(4) 


phase shift can be allowed for, so that cr can he assumed to 
have the form 

(T=e(a cos ktj +6 sin kj/) ( 8 ) 

Substitution of these sinusoidal relations into the general 
form of the boundary conditions, equations (5) , yields, after 
rearrangement. 


where m = UaIU. 

On the downstream side of the shock the velocity perturba- 
tions in the directions of J and ij, respectively, are (fig. 4) 

w„=(C/+dC7) cos (^+<r) + (F+dF) sin — 

> (5a) 

wj =~{U-\-dU) sin (^+<r) + (F+(fF) cos ) 

Equations (3) and (4) may be used to evaluate the right- 
hand side of equation (5a) . A first-order approximate result 
is obtained by taking cos c= 1 , sin <x= a and neglecting <r tan <p 
and <r cot in comparison with unity. It will be useful also 
to introduce the geometrical relation Ua = Wa cos 6, the 
definition UAlU=m, and to eliminate 8 by means of the 
oblique-shock relation tan (p=m tan 6. The final rearranged 
result is 


IF 


1—2 m + m^^ sin <p- 


7+1 


Ks( 


, sin^ <p~\ . 

cos ij5-| cos /C7J + 


COS ip_ 


— (l — 2 m + m^ 

_m \ -.11 


7+1 


w . 

Tu 


sin (pj sin ki? 

7 L 7+1 / cos <p 

a(m— 1) cos ip+ 2 J-m^sin ^sjcos kt; + 


b / 1 I 3 — 7 \ 

_ mV 7+1 / 


sin'= (p 
cos <p 


\-b(m — 1) cos (p sin Kij 




H9) 


V=-(wA-i ^+] 

tan ^+mo-^ sin tp 
tan ^+mcr^ cos <p — c sec ip 


w, 

u 


(5) 


These are the desired boundary conditions in a somewhat 
general form. 

In the present problem the perturbation Wa is associated 
with an incident sinusoidal shear wave parallel to Wa (or to 
^.i) (figs. 2 and 3). It will be shown later that a refracted 
sinusoidal shear wave parallel to W (or to f) will also arise. 
A suitable defining equation for Wa is 


— = i cos kriA (6) 

V\ A 

where k is the wave number (2x/A;=wave length). The 
corresponding argument for the refracted shear wave will 
involve ri and an altered wave number k. The arguments 
of the upstream and downstream waves must match along 
the shock front, so that 


Equations (9) give, in final form, the conditions imposed 
by the shock wave on the components parallel to ^ and rj, 
respectively, of the perturbation velocity immediately 
behind the shock; the parameters a and b therein governing 
the shock inclination tr are undetermined. These equations 
constitute the boundary conditions for the perturbation flow 
downstream of the shock. 

Evaluation of rotation term in governing equation. — Before 
equation ( 1 ) can be solved, the vorticity term (rotation 
term) on the right-hand side must be evaluated for the 
region behind the shock. A corresponding term has been 
evaluated in reference 7 for the flow behind a normal shock 
perturbed by an isoenergetic upstream disturbance. This 
work has proved a useful guide, but it has been necessary to 
make modifications both for the variation in energy (that is, 
in total enthalpy H) and for the inclination of the shock in 
the moving frame of reference. The derivation is as 


follows ; 

Downstream of the shock, the enthalpji H and the entropy 
s (and hence the vorticity) are constant along streamlines, 
and in the linear theory the streamlines are approximated by 


Thus, and ^ may be evaluated at 
’ dv bt] ■’ 


lines r) = constant. 

the shock and the result will hold downstream thereof if 
expressed as a function of r> alone (| eliminated) . 

The total enthalpy upstream and at the shock is 


kr)A~Kri along shock 
By geometry (fig. 3), Thus 

cos K 17 along shock ( 7 ) 

W A 

Since the disturbance is sinusoidal, the shock inclination <x 
can likewise be expected to be sinusoidal. For generality a 


H — Cj,Ta-\-~2 (l+A+W/i)^ 


Hence, at the shock 


bH 

bn 


shock 


( 10 ) 
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The entropy upstream of the shock is constant by virtue 
of the assumption of constant pressure and density there. 
The entropy change in crossing the shock is given in terms 
of the upstream velocity by (ref. 8, eq. 144) : 

In (Wa+WaY cos^ (« + <^)-^]x 

ny-l)(W^WAr cos^ (g + <r) + 2 ~[n 

L (7+1)(1^^a+Wa)^ cos^ (e + ff) J) 

Hence, on writing the differential and expanding the result 
under the assumption that Ta is constant and Wa/Wa and fl- 
are small, there is obtained 

Ss (m — ly a t&n 0^ I (Ha) 

and ^along shock 

Recall now that the governing equation (1) roads 


results will illustrate important features of the general case 
as well as provide a limiting case of the general solution, 
useful as a check. 

The horizontal wave is obtained by setting e=<p—0 in 
the earlier equations; as a consequence F-^0, W-^U, 
ij— > 2 /, and K->k. The governing equation reduces to 


;8Vxx+lAn/=— 1) sin h/ 


(14) 


where 


/?2=l-(/2 

The boimdary conditions (eqs. (9)) reduce to 


3 

Ue 


^=^=(m — l)u cos ky-\-{m—l)h sin kyy 
If Ue / 


(15) 


Particular integral and complementary function. — A par- 
ticular integral of equation (14) may be obtained by inspec- 
tion as 

(2m— 1) sin ky 


where the right-hand side is tlie rotation term in (piestion. 
The factors /?, and .s, have been evaluated in equations (10) 
and (11), respectively; substitution with use of the geometri- 
cal relations of figure 3 j'ields 


( 1 — + 


2 COS <p d 

cos^d i>ri 




U COS <fi{m — 1 ^ I®-!! ^ (12) 


where the right-hand side is to bo evaluated along the shock 
{x = 0) and expressed as a function of ?j alone. 

In the form (12) the governing equation has not yet been 
specialized to a shear flow that is sinusoidal. The substi- 
tution of equations (7) and (8) for WaIWa and <r, respectively, 
introduces the sinusoidal character; furthermore, the relation 
tan ^=m tan 6 can be used to eliminate 6; after simplifica- 
tion 

( 1 — 



fjYl 2 

sec ^-l-2(m— 1) cos <p-\-a - — ^ sin <fi 


'] 


Sin KT] 


, ^ (m— 1)=^ . 

+ ko ^ ^ sin (D 
m 


cos K7J 


X 

j (13) 


Equation (13) is the partial differential equation to be 
satisfied by the flow downstream of the shock subject to 
the boundary conditions (eqs. (9)). 

SOLUTION FOR HORIZONTAL WAVE 

The governing equation and boundary conditions have 
been set up for the general case of an inclined shear wave. 
It will be worthwhile to solve first, however, the much 
simnler special case of the horizontal shear wave. The 


'I'o obtain a complete solution there must be added a com- 
plementary function satisfjnng equation (14) with the right- 
hand side set equal to zero. The boundary conditions at 
x=0 require that the function possess a sinusoidal variation 
with y. Such a solution will also contain an exponential 
factor, showing either amplification or attenuation of the 
disturbance with distance x downstream of the shock; the 
case of amplification must be ruled out as physically unac- 
ceptable. These considerations limit the solution to the 
form 

4ic=Lfede ^ sin ky 

where d is a constant of integration. 

The complete solution is the sum of 4 'p fiud 

}j/= Ue ^ + sin ky (16) 

Evaluation of undetermined constants. — The velocity com- 
ponents are obtained from equation (16) as 

u = -4/y=Ueid^m—l-\-kde cos ky 
_kx 

v = — ^^\kx=^Uekde ^ sin ky 



The undetermined constants a, b, and d are evaluated by 
setting x=0 and comparing with the boundary conditions, 
equations (15), equating the respective coefficients of sin ky 
and cos ky. The.results are 


a = 0 


b 

d 


4|8m 
( 7 -|-l)(m-l) 


4m 


( 18 ) 


^(y+1) 
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u 

.. \ _ 

^ V 



?' h 


Z’ — 


i 


V 

1 Xz 

' , \ 


^ \ 
H U 


^ 

-A 

— •) 

•V K 

1 h 


— y 

~7 


t=!7f 2m — : 


4m 




alla;^0 (19) 


V — — ^JJf. — |— f e ^ sin ky 

7 + 1 


Just behind the shock 


Uo= Ue 


(-1 + 2 ^m) cos ky) 
4m 

U,^-,inky j 


4m ^ 

Vo= — pUf sin ky j 

and far downstream 

\ 

u„=Ut{2m— 1) cos ky) 

>x= 00 (19b) 

®„ = 0 ) 

These results and the associated streamline pattern are 
exhibited pictorially in figure 5. 

These perturbation velocity components downstream of 
the shock are to be compared with the corresponding 
velocity components in the shear flow upstream of the shock 
(cf. eq. (7)): 

Ua= Ua^ cos ky ) 

= U tm cos ky) 

Va=0 


The I’atio of u„!ua is 


^=2-— 
Ua m 


Since m=t+/f7>l in order that a normal shock exist, it 
appears from equation (20) that the normal shock always 
amplifies the horizontal shear wave, the maximum amplifica- 

tion of : being approached as the initial Mach number 
7+ 1 

approaches infinity. 


Shock perturbation. — The local inclination of the shock 
from the vertical is given by equation (8). With the previ- 
ously determined values of a and h (eq. (18)) inserted, and 
with ky in place of kij, the inclination is 


4|3em 

( 7 +l)(m- 


sin ky 


If the local shock displacement in the x-direction relative to 
the mean shock plane is called 8x(y), then 


Shock 

Fic.cre 5. — Pa.ssage of horizontal shear -wave through normal shock, 
sho-n'ing perturbation of shock and final amplification of shear 
wave. 

Velocity components. — Insertion of the value of d into 
equations (17) yields the final result for the velocity com- 
ponents valid everywhere downstream of the shock 


5a: = J* <r dy 


(21) 

Thus the shock displacement curve is in phase with the 
velocity perturbation in the shear wave upstream of the 
shock (fig. 5). 

SOLUTION WHEN FLOW DOWNSTREAM OF EQUIVALENT 
OBLIQUE SHOCK IS SUBSONIC (W<1) 

The present case is a generalization from the horizontal 
wave just discussed to a wave of arbitrary inclination 6. 
The restriction to a subsonic mean velocity W behind the 
equivalent oblique shock insures a qualitative similarity of 
the flow: the governing equation is elliptic in both cases. 
Accordingly, the horizontal-wave result can serve as a guide. 

Governing equation and particular integral. — The govern- 
ing differential equation (13) may be written in abbreviated 
form as 

— nUt{A sin ki]—B cos kij) (22) 

where 

A=r sec (p + 2(m — 1) cos <p+a gin y, ^ 

L m J# 

I (23) 

m I 

A particular integral is seen to bo 




A . B 
— sm KTi cos Kt) 

K K y 


Complementary function.— From the result for the hori- 
zontal wave the complementary function should be expected 
to attenuate exponentially downstream of the shock, and 
from physical considerations the attenuation should depend 
upon the distance measured normal to the shock front, that 
is, upon X rather than, say, ^ The functional form that has 
the desired attenuation and possesses a sinusoidal behavior 
at the shock is 


c. —a^io COS ip —rj sin V) 


or [a(f sin r/ cos <p)] > (25) 

cos ) 


where | cos (p—rj sin (p may be recognized as just 


I 
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The arbitrary constant a in equation (25) is determined 
by a consideration of the boundary conditions (eq. (9)): the 
argument of the cosine must reduce to kt? along the shock 
front, where |=ij tan <p. This requirement gives a=K cos <pl0‘. 
Finally, when constants of integration c' and d' are included, 
the complementary function is written as 

, -tt -^oos*>({cos*>-:isin*>) (r K COS . , 

\l/a=Uee ^ <1 c' cos — sin^jd- 

PJrt cos sin ^ sin <p+^Jrt cos ^)J j (26) 

Velocity components with undetermined constants. — The 


complete solution for the perturbation stream function is 

^='/'p+'Ac 

This expression (cf. eqs. (24) and (26)) contains four arbitrary 
parameters a and b (which occur in A and B, respectively, 
eqs. (23)) and c' and d', which remain to be determined. 
First the corresponding expressions for the velocity compo- 
nents will be obtained — they will be needed anyway — and 
then the boundary conditions on these velocities at the shock 
wave will be applied for the determination of a, b, c', and d'. 

The perturbation velocity components in the direction of ^ 
and V are and w'= — | 8 „Vf, respectively; by differentia- 

tion of equations (24) and (26) there results 


W . I 73 • I ^ cos *> ({ cos V-I, sin v) 

cos sin Kv+P e ^ 

U e 


, - I j • \ • K cos ip (I sin + cos v?) 

(— cos <p+d sm <p) sm — — 


W' . cos (f cos v—n sin (»); 


|^(c sin (p-{-d^a cos (p) cos 

] 


K COS tp ($ sin (pA^Jn cos tp) , ') 


t / a 1 ja • \ K COS <p sin <p A ri COS <p) 

cos (p—dPw sm <p) cos — 


, _ . , Jo 2 \ • K COS ip (J sin ?) cos ^)' 

(c/3„ sin p + d/Sn. cos <p) sm 


(27) 


where c' and d' have been absorbed for convenience into 
new constants 


c = c'K^m cos tp 


d = dAfia cos <p 


CEADF^ 
® ^ C'^AD^ 

, CF-DE 
b-m 

c,=±D'-F' 

m 


(29) 


The undertermined constants may now be considered as a, b, 
c, and d. 

Conditions along the shock on the downstream side have 
been designated by the subscript zero; here $ cos p=j? sin ip, 
and the arguments of the exponential and sine and cosine 
terms reduce to zero and kt?, respectively: 

sin <pA~^ cos tp^ cos 107 -f- "| 


where 


d=- D' 
m 


tan p-[(m-l)^-b 1 


2 (m — 1 )' 

7 + 




sm <p cos (f 


— l+(m — l)cosV D' 


(30) 


/3„ cos <pA-p sin sin kjj 

cos tp—dfiy, sin <p) cos K-qA 
/3~^(c/3„ sin ipAdpJ^ cos A) sin ktj 


(28) 


Evaluation of undetermined constants. — Equations (28) 
must agree identically with the boundary conditions (eqs. 
(9)) imposed by the shock wave on Wo and w„' . Therefore 
the respective coefficients sin Kq and cos Kq are to be equated; 
this yields four simultaneous equations for the four undeter- 
mined constants a, b, c, and d. In the reduction of the solu- 
tion to final form certain alternative forms of the oblique- 
shock relations, given in appendix C, have been used. The 
results are 


»i)+2(m-l) 

1 ) sin V cos pj=|l P' 


SOLUTION WHEN FLOW DOWNSTREAM OF EQUIVALENT 
OBLIQUE SHOCK IS SUPERSONIC (W>1) 

When the mean velocity W behind the equivalent oblique 
shock is supersonic, the solution must exhibit Mach waves. 
If the cross-stream velocity V of the moving reference frame 
is subtracted out, these waves appear to be moving down- 
ward (cross-stream) with the velocity V. If another trans- 
formation of axes is made so that the reference frame is 
“convected” downstream with the stream velocity U, then 
the Mach waves can be identified as plane sound waves 
moving normal to the wave fronts with sonic velocity. Mach 
waves and plane sound waves are, of course, the same 
phenomena viewed relative to different frames of reference. 


OQaRJO— Fi4, 



8 


REPORT 1164 NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Governing equation and particular integral. — The govern- 
ing equation ( 22 ) changes from elliptic to hyperbolic when 
W exceeds unity (that is, when W is supersonic). The 
particular integral is unchanged thereby and is still given 
by equation (24). It is found that the final solution yields 
b=0 (and hence B=0), and so it is convenient to delete the 
S-term at the outset; the particular integral is thus 


are w=^i} and respectively; the expressions are 


w 

Ut 


-[a cos 


c sec (f> k(^+Pv,v) 

C/OS 


-f-tan Ip /3u,-ftan 




( 32 ) 


w' |3ttC sec <p k(^-|-/ 3 m 7 ?) 

Ue /3„-f-tan jS„-l-tan ^ 




UiA 

K 


sin Kt] 


Complementary function. — The complementary function 
satisfying equation ( 22 ) must be of the general form 

where The function/ represents Mach waves 

inclined downward by the Mach angle p from the ^-axis and 
the function g represents Mach waves inclined upward by the 
Mach angle. If attention is restricted to the range of shear- 

wave inclinations 0<B<~! then the p-family of Mach waves 

/ii 

can be shown to represent disturbances overtaking the shock 
wave from behind. This property is related to the fact that, 
for a finite shock strength, the Mach angle is always greater 
than the angle between the shock and the ^-axis. Since the 
disturbances actually originate at the shock wave by virtue 
of the passage therethrough of the initial shear wave, such 
Mach waves cannot arise, and the p-function must be zero. 
In what follows it will suffice to limit the discussion to the 

specified range since the results for the remaining 

range O<0< — ^ are readily obtained therefrom from sym- 

me try considerations. 

The function / must reduce to 

/ ~sin Ki\ 

along the shock front, where ^=?j tan <p, in order to satisfy 
the boundary conditions (with 6=0). A suitable comple- 
mentary function is therefore 


i>c= 


Vec" 

K 


sin 


;3„-|-tan <p 


where c" is a constant of integration. 

The complete solution for the perturbation stream function 
is thus 

-A=vJ'c+-/'c=^[Asin «, + c"sin-g^^] (31) 

This expression contains two arbitrary parameters a (occur- 
ring in A) and c" which remain to be determined. First the 
corresponding expressions for the velocity components will 
be obtained, and then the boundary conditions on these 
velocities at the shock wave will be applied for the deter- 
mination of a and c”. 

Velocity components with undetermined constants. — The 

perturbation velocity components in the direction of ^ and 17 


where the constant c” has been absorbed into a new constant 
c=c" /3„ cos <p. The undetermined constants are now a and 


Along the shock ^=17 tan <p, and the arguments of all cosine 
terms reduce to kti; the expressions for the velocity com- 
ponents w and w' become 


Wq 

Ui 

Wo 

Ut 


G 


c sec p 
Pw “k tan p , 


3wC sec p 
j8„-ktan p 


cos fC?7 




(33) 


Evaluation of undetermined constants. — Equations (33) 
must agree identically with the boundary conditions (eqs. 
(9)) imposed by the shock wave on Wo and w„' . If the respec- 
tive coefficients of sin K 17 and cos kt; are equated, there results 
6=0 and two simultaneous equations for a and c. Thus, the 
initial specification of 6=0 has been justified a posteriori. 

The solutions may be written in the form 


_ G' + GF' 

c=^D'-F' 

m 

where 


(34) 


6" = 2^^4m-2[l-k(m— 1) cos' p] 
7-f-l 

D' ={m — 1)[1 -|-(m — 1 ) cos' p] 


(35) 


F' = 2(tn — 1 ) sin p cos p 
g^i ^g/^^ =tan (p-p) 

; 8 „-|-tan p 


where /.t = cot '/3„, is the Mach angle. (The definitions for D' 
and F' herein are unchanged from those included in eqs. 
(30).) 


RESULTS AND DISCUSSION 


VELOCITY FIELD 

The velocity field downstream of the shock wave, produced 
by convection of an oblique sinusoidal shear wave through 
the shock, has been calculated; the results are distributed 
through the preceding sections. The main results will now 
be presented in more compact form, simplified to aid in the 
geometrical interpretation. (The special case of the hori- 
zontal shear wave was discussed earlier.) 

Frames of reference. — The analysis has been carried out 
in a special frame of reference in which the flow is steady; all 
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Figure 6. — Relative motion of reference frame moving with general 
downstream flow (convected frame) and reference frame of analysis 
(steady-flow frame). The steady-flow frame moves downward 
along the shock front with a component velocity V and carries the 
ripple pattern with it. 

formulas will be given relative to this steady-flow frame. 
Also of considerable interest is a frame of reference convected 
by the mean flow downstream of the shock; this frame is at 
rest relative to the general mass of fluid there. The relation 
between the two frames is shown in figure 6. Formulas rela- 
tive to the steady-flow frame may be converted to apply to 
the convected frame by means of the transformations 


v^v 

> 

Ut 

y^y+Vt. 


(36) 


The criterion on W. — Although the stream velocity U 
downstream of the specified normal shock (fig. 1) is always 
subsonic, the nature of the flow depends primarily on the 
stream velocity W downstream of the equivalent oblique 
shock (figs. 2 and 3), which may be either subsonic or 
supersonic. The velocity W may also be interpreted as the 
relative velocity of the steady-flow frame of reference and 
the convected frame (fig. 6) . Two forms of the solution for 
all flow quantities thus appear, one for the su^onic range 
W<1, the other for the supersonic range VF>1. The 
dividing line W= 1 is what has been designated “the criterion 
on W” at the head of this section. Since W depends on the 
initial Mach number Ua. and the inclination 6, the equation 
W=1 gives, in effect, a relation between a critical value of 
d and TJa- The relation is conveniently expressed in terms 
of m=UAlU, which depends on Ua (see appendix C) : 


Ocr= ±tan 


V 






(37) 



s-l 


CD 

C 

o 


II 

6 


Figure 7. — Upper curve shows variation with initial Mach number 
of critical wave inclination for which W is sonic. Lower curve shows 
variation of in= (7^/17= tan ^/tan 9 with initial Mach number. 


A graph of \dcr\ against Ua is given in figure 7. 

The equations that follow, both for W<1 and for W>1, 

are in a form that applies for O<0< ^ ^whence 0<<p<^^- 

Corresponding relations for— ^<^<0 can be 

obtained by use of the symmetry of the flow with respect to d 


<p- 


Resultant velocity, W<\. — Equations (27) may be recast 
in the form 


\wa\ 


=S cos [k„(i/ — X tan 95)-f-6s] + 


n(x) cos [K„(t/ — X tan ^') + 5i,]V 


\ 


= /3„n(x) sin [K,{y — x tan <p') + 5p] 




w 

where ® 

= W)i6= amplitude of sinusoidal velocity Wa in initial 
shear wave 
Ky—K cos tp=k cos 0 


A = A{a), B=B{b) 

Tft 

' ^ m iS 
5j=tan-i 


* Typographical errors appearing in Technical Note 2864 (which is superseded by the 
present report), namely, in Sp (W<1), and in n and <p' (W>l), have been corrected. These 
errors were not embodied in the numerical calculations and curves of the present report. 

® In all the arc tangent designations, the quadrant of the angle is to be determined by the 
respective signs of numerator and denominator; for exam ple, 5«= tan~K~-S/y4. ) is to b e re- 
garded as determined by the joint conditions sin — cos 8» =AiyA^+B>. 




10 


REPORT 1164 — MATTON AT, ADVISORY COMMITTEE FOR AERONAUTICS 


(p'=—ian~ 


, tan <p 


Sy=tan ^ 


c^g—d tan <p 
dPg-\-c tan (p 


The functions A and B are given by equations (23) and a, b, 
c, and d are the initially undetermined constants which have 
been evaluated in equations (29) and (30). 

Resultant velocity, W>\. — Equations (32) may be recast 
as follows : 


U) 

cos Ky{y—x tan <ff)+n cos Ky{y—x tan 
— | = j8«, n cos Ky(y—x tan <p') 



where 

= lTA€=amplitude of sinusoidal velocity Wa in initial 
shear wave 




cos 6 
m 


A; 


A = A(a) 


Ky—K cos <p=k cos 8 

„ cos 8 sin M 
II — c f \ 
m cos(^— ju) 

<p'=(p—fi 

M=Mach angle=cot“* /3a, 

The function A is stUl given by equations (23), and a and c 
arc evaluated in equations (34) and (35). 

Shear-wave component. — The cosine in the <S-term is con- 
stant along lines y—x cot ¥5=constant; such lines are inclined 
at an angle (p with the horizontal and are thus parallel to 
the |-axis. Since w is parallel to J and w' is parallel to v, it 
is seen that the S-term represents a pure shear flow parallel 
to the |-axis. Stated otherwise, this is a rotational flow; 
the rotation (or vorticity) is just S2, which was evaluated 
earlier in terms of gradients of entropy and total enthalpy 
(cf. eqs. (1) and (13)). The shear flow may be described 
also as an incompressible, plane, transverse, sinusoidal wave. 

The amplitude and phase of the shear wave are compared 
with those of the initial shear wave in figure 8 for an initial 
Mach number of 1.5. The amplitude amplification ratio 
is S and the angle of phase lead is ; both are plotted against 
the initial wave inclination 8. Th_ere is seen to be a small 
phase lead in the subsonic range (1T<1) and none at all in 
the supersonic range (1T>1). The amplification is nowhere 
less than unity, with a cusp-like peak of 1.73 at the sonic 
point W—\. 

Pressure -wave component. — The remaining terms in equa- 
tions (38a) and (38b), involving the factor II, correspond to 
an irrotational velocity field, or potential flow. That is, if 



Figure 8. — Amplification and phase shift of velocity in shear wave on 
passage through shock. Initial Mach number Uas 

the derivation is traced backward, the Il-terms are found to 
have come from the complementary function, which is a 
solution of the governing equation with the vorticity 0 set 
equal to zero. This part of equations (38a) and (38b) 
defines what may be called a pressure wave since there is 
associated with it a first-order pressure field : the shear wave 
contributes nothing to the pressure. 

The pressure wave may be interpreted as a distribution of 
sound waves. This interpretation is particularly evident 
for the case TTj>l, where the solution has been obtained in 
the form of Mach waves: if a transformation is made from 
the present special frame of reference, relative to which the 
flow is steady in time, to a frame moving with the general 
stream, then the Mach waves will reappear as plane sound 
waves moving normal to themselves with sonic speed. 

The same transformation results in somewhat more com- 
plication when TT<(1: the resultant pressure pattern does 
not then propagate with the speed of sound, but it can be 
represented (as can any two-dimentional irrotational gas- 
flow field) as a superposition of cylindrical sound waves 
which individually propagate with sonic speed. The associ- 
ated velocity pattern in this case exhibits the following 
features, which are brought out by an examination of equa- 
tions (38a) : The radius vector in a graph of w' against w 
(hodograph) moves in an ellipse when x is held fixed and y 
varied; the major and minor axes are YI\wa\ and /3„,n|w^|, 
respectively. At a:=0 the phase angle relative to the inci- 
dent shear wave is 8p. On the other hand the argu;ment of 
the cosine and sine is constant along lines y—x tan 8'= con- 
stant; these are lines inclined at an angle 8' to the horizontal. 
Along such lines the perturbation velocity {w, w') re.mains 
constant in direction but attenuates exponentially with x; 

the exponent is — -2 • 
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FifjuBE 9. — Amplitude and phase of velocity components in pressure 
wave generated_by passage of shear wave through shock. Initial 
Mach number Va, 15- Parallel line.s in inset figure are lines of 
constant pha.se. 

For the case 1F>1, the velocity pattern associated with 
the pressure wave is much simpler (eqs. (88b)). The per- 
turbation velocity vector {w, w') is constant along lines 
y-x tan ^' = constant and is, in fact, normal to such lines. 
In this case = Mach angle, and these arc just the 
Mach lines (or envelopes of the sound waves); they are 
inclined downward by the Mach angle p relative to the 
J=axis. It will be noted that the definitions of tp' , the 
inclination angle of the lines of constant phase, agree at 
W=l, although expressed differently for W<C 1 and for TF>1. 

The amplitude and phase of the w and w' components of 
the velocity in the pressure wave are compared with the 
amplitude of the initial shear wave in figure 9 for an initial 
Mach number of 1.5. The amplitude amplification ratios 
are II and /3„n, respectively; n, /3„n, and a phase angle (lead) 
Sj, are plotted in the curves against the initial wave inclina- 
tion d. In the subsonic range (PF<l) II and /3„n attenuate 
exponentially with x and only the values for a;=0 are plotted. 
The phase lead varies from 180° to zero in this subsonic 
range and remains zero throughout the supersonic range 
(1F>1). A rather striking feature is the relatively small 
perturbation velocity in the supersonic range. Thus, 
although the incident shear wave can give rise to a simple 
sound wave upon passing into the shock wave, the particle 
velocity in this sound wave amounts for most cases to 10 
percent or less of the velocity in the initial shear wave for 
Ua^I.5. 

PRESSURE FIELD 


Here Wp is that component of the perturbation velocity 
associated with the pressure wave and directed parallel to 
W (that is, along the f-axis). Equation (Bll') may be 
recognized as the linearized Bernoulli equation as limited to 
the velocity in the pressure wave. 

Upon substituting for W and W and using for Wp equations 
(38a) and (38b) with the <S-ter.ms omitted, there results 


V Va 


2y7nB. sec <p 
( 7 -f l)m— ( 7 — 1 ) 


cos [Ky{y — x tan <p')-\~Sp\ 


(39) 


where Sp is to be taken as zero in the supersonic range of W. 
This result for the perturbation pressure is proportional to 
n sec <p] n has been plotted in figure 9, together with 6p, as a 
function of wave inclination d for Wa = 1-5. 

DENSITY FIELD 

The density perturbation is related to the velocity and 
entropy perturbations according to equation (B12) of appen- 
dix B; in the present notation this is 


(B120 

P W Cp ^ 

The term in Wp is the contribution of the pressure wave. 
This term differs fro.m 6p/p (eq. (Bll')) by a simple factor I/7, 
so that the contribution is obtained at once fro.m equation 
(39). 

The term in Ss is the contribution of the shear wave. The 
entropy perturbation Ss has not been given explicitly before, 
but it can be obtained from equation (11a) by use of geomet- 
rical relations and the known result for o- (see following sec- 
tion). Upon evaluation, the term in 5s is found to be 


\Wa\ 
^ p V A 


2{m—\y cos B. a I 

- — ,-^j [{a tan 6—1) cos 

T+ 1 , 

- — m— 1 
7— 1 


h tan 6 sin Ky] (40) 


TEMPERATURE FIELD 

The temperature perturbation follows at once from the 
pressure and density perturbations according to the equation 
of state; the appropriate small perturbation form of the 
equation is 

T p p 

SHOCK-WAVE PERTURBATION 

The local perturbation in the shock inclination angle may 
be written (cf. eq. (8)) 

<T=t(a cos Kyy-\-b sin Kyy) 


It is shown in appendix B that the perturbation pressure is 
related to the velocity according to equation (Bll); in the 
present notation this becomes 

dp= — pWWp \ 


where a and b are evaluated in equations (29) for TF<(1 and 
equations (34) for W)>1 (6 = 0 for W)>\). 

The local shock deflection 5a: from the plane x=0 is obtained 
by integration of the slope o-; 


P 


w 



or 


(Bll') 
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Figure 10. — Amplitude and phase of ripples developed in shock by 
passage of shear wave. Initial Mach number Ua, 1.5. 


The result may be put in the form 


5x = \wa\ 


X Va^-f 6^ 
2irUA 


cos (Hyt/ “|- Sghock) 


(41) 


where 5shock=tan”‘ phase angle and X=2ir//: is 

the wave length of the initial shear wave. 


For a given wave length the factor is proportional 

to the amplitude of this sinusoidal corrugation in the shock 
wave; is plotted against the initial wave inclination d 

in figure 10. The phase angle 4h»ok is also plotted: the 
shock-wave corrugation is in phase with the initial shear 
wave (5,h„„k=0) when the initial wave is horizontal (6=0). 
The shock corrugation progressively lags the initM shear 
wave as 6 is increased until the sonic condition 1F=1 is 
reached; at this point the lag is 90°, and this value is main- 
tained throughout the range 1F>1 as 6 is increased to 90° 
(vertical initial shear wave). At 0=90° the amplitude 
factor has fallen to zero: a vertical sinusoidal shear 

wave passing by convection through a vertical shock wave 
causes no perturbation of the shock form or position. 


INTENSITY or SOUND FIELD 

The analysis implies that the interaction of turbulence with 
a shock wave does not give rise to any great amplification of 
the fluctuation energy, but it does provide a transformation 
from a relatively quiet form (initial turbulence) to a relatively 
noisy form (final turbulence plus sound field). On an 
acoustic scale the level of the noise generated is found to be 


relatively high.® It will sufiice for an order-of-magnitude 
estimate to replace the turbulent field by a single plane wave, 
or Fourier component, with the same kinetic-energy density. 
Roughly this implies that the root-mean-square turbulent 
velocity is to be identified with 0.707 |wa|. 

The sound pressure is proportional to n sec <p, where n is 
plotted in figm’e 9. The relatively high values indicated 
for the subsonic range attenuate rapidly with distance x 
downstream of the shock; when x appreciably exceeds several 
wave lengths, the values are negligible compared with those 
in the supersonic range. A rough average over all wave 
inclinations 0, assuming the subsonic range contributes 
nothing, gives |II sec ip\ ^0.082; this value will be used in the 
noise estimate. 

The noise level in decibels relative to the standard reference 
level 6p„=2X10"'® atmospheres is given by 

db — 20 log ^ 

^ Spo 

where the 8p’s are root-mean-square values. By equation 
(39) 

5p rms \wa\ 

~i> in 

As an example the noise level generated by the turbulence 
passing through a normal shock in a representative super- 
sonic wind tunnel will be estimated. A root-mean-square 
velocity of turbulence of 1 foot per second is assumed to exist 
in the test section where the mean speed is 1400 feet per 
second and the Mach number is 1.5 (174=1.5). Thus 
0.707 |w 4| and Ua are taken to be 1 and 1400 feet per second, 
respectively. A summary of these and the remaining para- 
meters of the example is 


V2 ym 


l)m — (7— 1) 


n sec tp 


(43) 


I n sec ip I =0.082 
0.707 |w 4 | = l foot per second 
(74=1400 feet per second 
m = 1.862 (~04=1.5) 

7=1.4 

p = 0.670 atmosphere (~1 atm. reservoir pressure) 
5p„=2X10“‘® atmospheres 

The estimate based on equation (43) gives a pressure 
perturbation 8p rms/p=7. 50X10“®, and by equation (42) 
the corresponding sound pressure level is 108 decibels. This 
represents very intense noise, reaching a level which can 
damage the ear on continued exposure (ref. 9). This noise 

« This is more a testimony to the great sensitivity of the ear than to the energy content of the 
noise; thus the ear-splitting noise in the vicinity of a jet engine (say, 140 decibels) is produced 
by pressure fluctuations of ±0.003 of an atmosphere. 
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Figure 11. — Construction for translation V to render either of two 
sound-wave patterns stationary in a main stream Ua- 

estimate is thought to be conservative, corresponding to a 
supersonic wind tunnel with a relatively low level of turbu- 
lence. It appears probable that many tunnels will consid- 
erably exceed this level. 

GENEHALIZATION TO OBLIQUE SHOCKS 

The analysis refers to flow through a normal shock, but the 
results are easily generalized for oblique shocks. In the 
oblique-shock case the component of the upstream velocity 
normal to the shock plaj^s the role of Ua ; the component par- 
allel to the shock is ignored in formulating the equivalent 
steady-flow problem. A formal approach is to retain the 
present definitions wherein Ua is the actual upstream veloc- 
ity (taken horizontal) and 6 and <p are referred to the horizon- 
tal; the oblique shock is assumed inclined by some angle a 
measured clockwise from the vertical. Then the present 
formulas will be generalized to apply to the oblique shock if 
the following transformations are made: 

Ua—^Ua cos a 
6 — ^6 -j- a 
a 

RELATED PROBLEMS 

The sound field produced downstream by the convection 
of turbulence through a shock has been discussed. Also of 
interest are sound fields incident upon a shock in the absence 
of turbulence. The elementary sound disturbance is the 
plane sinusoidal wave: a longitudinal wave. The passage of 
such a wave through a shock, which is an unsteady-flow 
problem, can again be converted to an equivalent steady- 


flow problem by transformation to a reference frame moving 
with a suitable velocity parallel to the shock front; in this 
frame the sound-wave pattern will appear as a stationary 
Mach wave pattern. A diagrammatic construction is shown 
in figure 11. Note that either of two sound patterns of 
uniquely related inclinations may be rendered stationary by 
a given choice of V; the two patterns may be identified with 
the two families of Mach waves in a stream of supersonic 
velocity Wa- 

The equations for the boundary conditions at the shock 
and the vorticity behind the shock will be modified from 
those for the present case of the shear wave, but the general 
character of the solution will be unchanged. Thus, a shear 
wave as well as a sound wave will appear downstream of the 
shock. The discussion will be carried no further here; the 
solution has been obtained in reference 5 by the unsteady- 
flow method. 

The interaction of a sinusoidal Mach wave pattern with a 
normal shock constitutes a simple special case: here the 
velocity V of the moving reference frame may be taken to 
be zero. This problem lias been solved in general terms by 
Adams (ref. 7); he limited his discussion, however, to the 
vicinity immediately downstream of the shock. The char- 
acter of the flow further downstream can be inferred from the 
parallel that exists between this problem and the problem 
herein of the horizontal shear wave: in both cases V is zero. 
The asymptotic flow far downstream is therefore a horizontal 
sinusoidal shear wave. Near the shock the wave is modified 
by transveme and axial components (with associated pres- 
sure perturbations) which attenuate exponentially with dis- 
tance downstream of the shock (cf. fig. 5). 

According to these considerations, sinusoidal corrugations 
in a wind-tunnel wall, or a plate, upstream of a plane shock 
wave will general, e a horizontal sinusoidal shear flow. Such 
a shear flow might have applications in special experimental 
work. 

CONCLUDING REMARKS 

The effects produced by the convection of an inclined 
plane sinusoidal shear wave through a normal shock have 
been analyzed. Such a wave may be interpreted as a single 
spectrum component of a turbulent field; that is, the turbu- 
lent field can be represented as a superposition of such shear 
waves of all orientations and wave lengths (Fourier integral). 

When the turbulence is convected through a shock, the 
individual waves do not mutually interfere if, as specified 
herein, the intensity is sufficiently low: thus the modified 
field downstream of the shock can be obtained in principle 
by superposition of the modified individual waves. In prac- 
tice the initial wave distribution is known only statistically, 
and statistical changes only can be calculated. In either 
case the present analysis of the behavior of a representative 
individual wave constitutes a prerequisite to the determina- 
tion of the changes in the weak turbulent field. 

It is found that a sinusoidal shear wave of arbitrary in- 
clination as it passes into the shock gives rise downstream 
to a shear wave of altered inclination and altered amplitude. 
In addition, there is generated a “pressure wave”: an addi- 
tional velocity field with associated pressure disturbances 
that can be recognized as sound waves. 


14 


REPORT 1164 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


The analysis is made in a frame of reference moving with 
a certain velocity W referred to axes at rest relative to the 
general mass of fluid downstream of the shock: W is the 
vector sum of the reversed downstream velocity and the 
cross-stream speed of the ripple pattern in the shock wave. 
The results depend crucially on whether W is subsonic or 
supersonic: when W is subsonic both the shear wave and 
pressure wave are shifted in phase relative to the initial 
shear wave, and the pressure wave shows an exponential at- 
tenuation downstream of the shock; when W is supersonic 
there are no phase shifts, and the pressure wave takes the 
form of a plane, undamped, sinusoidal sound wave. 

A weak initial shear wave is found to produce a surpris- 


ingly intense pressure wave or sound field downstream of the 
shock, as measured in acoustic terms. This implies that the 
convection of relatively low-intensity turbulence through a 
shock win generate a very intense noise field in the down- 
stream region. In an example the noise pressure level gen- 
erated by turbulence in a representative supersonic wind 
tunnel was estimated to be of the order of 108 decibels. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 26, 1952 


APPENDIX A 


SYMBOLS 


The 


A 

a 

B 

h 

C 

C 


c 

Cp 

D 

D' 

d 

d' 

E 

E' 

F 

F' 

6 

H 

k 

M 

m 

'P 

S 

s 

T 

t 

U,V 

U 

u,v 

W 

W 

w.w' 


parameter 


following symbols are used in this report: 
function defined in eqs. (23) 

parameter in shock-wave perturbation (eq. (8)); 

also speed of sound 
function defined in eqs. (23) 
parameter in shock-wave perturbation (eq. (8)) 
function defined in eqs. (30) 

function defined in eqs. (35) 

/ — c"j3«, cos ip for ^>1^1 
\ = c'Kj3„ cos <p for IPOl/ 
constant of integration 
constant of integration 
specific heat at constant pressure 
function defined in eqs. (30) 
function defined in eqs. (35) 
parameter ( = d'K/3w cos tp) 
constant of integration 
function defined in eqs. (30) 
function defined in eqs. (35) 
function defined in eqs. (30) 
function defined in eqs. (35) 
function defined in eqs. (35) 
stagnation enthalpy (per unit mass) 
wave number of shear wave in region A (incident 
shear wave) 

Mach number ((7/a, appendix B) 
velocity ratio across normal shock (UaIU) 
pressure 

relative amplitude of refracted shear wave (see 
eqs. (38a) and (38b)) 
entropy (per unit mass) 
temperature (absolute) 
time 

stream velocity components in x- and y-directions 
(fig. 3) (equivalent steady-flow problem) 

Mach number associated with U (Ula) 
perturbation velocity components in x- and y- 
directions, respectively (fig. 3) 
stream velocity in ^-direction (resultant of U and 
V) (equivalent steady-flow problem) 

Mach number associated with W (IT/a) 
perturbation velocity components in and jj- 
directions, respectively (fig. 3) 


x,y 

X 

b, 

S, 

f 

e 

K 

n 

p 

O' 

<p 

v' 

i' 

’Ac 


that part of w associated with pressure wave 
rectangular coordinates (fig. 3) 

■y/l-U" 

W<1 

\^W-i W>i 

ratio of specific heats 

phase lead of pressure wave relative to incident 
shear wave 

phase load of refracted shear wave relative to 
incident shear wave 

measure of strength of incident shear wave 
(M/Wa) 

inclination of lines of constant phase in incident 

shear wave (figs. 1, 2, 3, and 4) 

critical value of 6 for which W=1 (function of Ua) 

wave number of refracted shear wave 

Mach angle associated with W (sin”‘ {IjW)) 
inclined rectangular coordinates (fig. 3) 
relative amplitude of velocity component w in 
pressure wave (see eqs. (38a) and (38b)) 
fluid density 

perturbation in local shock angle (fig. 4) 
inclination of lines of constant phase in refracted 
shear wave 

inclination of lines of constant phase in pressure 
wave 

perturbation stream function 
complementary function (component of ^) 
particular integral (component of ^) 
vorticity (rj— Rj,) 


Subscripts: 

A region A (upstream of shock) 

o evaluated at shock, on downstream side 

X, y, L V indicate the corresponding partial derivatives 

^e. g., exception is Ky=K cos (p= 

k cos 6 


(Unsubscripted velocity components, pressure, and 
density refer to region downstream of shock.) 

Prolix 

6( ) increment in ( ) 
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APPENDIX B 


LINEAR PERTURBATION THEORY FOR ROTATIONAL FLOW 


The generalized governing equation for the stream func- 
tion can be obtained b}^ extending Sears’ constant-energy 
development (ref. 6) to include the effects of variation of 
energy (total enthalpy). A different approach is, however^ 
employed herein. Equations for the pressure and density 
fields are also obtained. 

In applying the results of this appendix to the develop- 
ments in the main text it is to be noted that the x- and 
?y-axes herein will go over, respectively, into the and 
ij-axes therein; this is a consequence of the difference in 
direction of the main stream in the two cases. There is a 
coi'responding change in the notation for the velocity com- 
ponents. 

GOVERNING EQUATION 

Basic equations. — Consider the steady two-dimensional 
adiabatic flow of an in viscid fluid \vith local velocitj' -a', 
])ressure p, densitj” p, temperature T, and entropy s. As- 
sume only small perturbations from a uniform horizontal 
flow such that u'=U+u, v'=v, with ulU, and 

also Spip, ipjp, etc. <(<(1. Then the basic flow equations 
may be linearized by neglecting quantities of order ujU, 
and so forth, in comparison with unity. A convenient form 
of these linearized equations is 

Continuitv: 'Ux+ rj,-}-- ^^ = 0 (Bl) 

■ p ut 


State: 


&p_ 


pa- 




E 


nergy: 


D.S 

Dt 


= 0 


(B2) 

(B3) 


— P jr=- pU UX~\ 

Momentum: > (B4) 

— Py—pUvx) 

where DjDt signifies the Lagrangian operator for differentia- 
tion following the fluid motion. 

Elimination of density from continuity equation. — The 
Lagrangian form of the state equation is, by virtue of the 
energy equation, 

1 

p Dt pa- Dt 


W -\-v)px-^vpy\ 
pO, 

Upon linearizing, assuming py and px to be of comparable 
magnitude, this is 

1 _ J-‘Vx 

p Dt pa? 


and by use of tbe first momentum equation 
1 Dp 


Dt 


= -MHx 


The linearized continuity equation (Bl) ma}' accordingly be 
written 

(1— r„=0 (Bo) 

Formulation pf governing equation. — Define a stream func- 
tion ^ such that 

u=\py 

v=^-{l-M-)i^x (B6) 

Then equation (B5) is identically satisfied by u and r as 
defined in equation (B6). The governing ecpiation for xj/ is 
now obtained by expressing the vorticitj' Vx—Uys9. in terms 
of xj/: 


( 1 — d /-) xlxxx -\-^yy=—9. (B7) 


A useful expression of the vorticity in terms of gradients of 
entropy and total enthalpj' is given in reference 10, equation 
(8.3), as 



2)H j, ds\ 
dp. bn) 


Avhere g is the resultant velocity and bjbn indicates differ- 
entiation normal to a streamline. In the small-perturbation 
flow the. streamlines are approximated bv* the lines y= con- 
stant, so that d/&7jsd/dy; also g^^V. Thus 


_o 



fbH rp d.s \ 
\by by) 


The governing equation for xj/, equation (B7), can now be 
amplified to read 

TT 

(1 - M^)^„-h Xl>yy= - fi 


This equation and its companion 




(B6) 


v==-{l-M-)xlxx 


constitute the simplified generalization of Seai-s’ governing 
equation for linearized rotational flow (ref. 6, eqs. (12) and 
(15)); Sears’ equation is restricted to flows of constant total 
enthalpy H. 

Equation (B8) exhibits the following veiy interesting 
properly: In the small-perturbation velocity field considered here 
the ejffect of the rotation or vorticity Q is independent of how it 
ai-ises, whether from a gradient of entropy or a gradient of total 
enthalpy, or a combination of both. And it is only through their 
contribution to {and perhaps to the boundary conditions) that 
variations in entropy and total enthalpy affect the velocity held 
at all. 
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PKESSURE FIELD 

Equations lelating tlie pressure distribution to the velocity 
distribution will now be derived: The momentum equations 
(B4) maj' be rewritten in the form: 


•> 

p 

^+Uu,^-U9. 

P J 

since Q^r^—Uy. 


(B9) 


Now let consideration be limited to special types of flow such 
that 


U = w'd- 


r=r 



(BIO) 


where u' , v' is an irrotational flow (r'^— a'j,=0) and 

u" = u"(y) is a pure shear flow parallel to r (r" = 0).' Then 

the vorticity 0 is given by 


and 






Thus equations (B9) become 


^+ Uu '^-=0 

P 


^+ Uu',^0 

P 

These two equations are equivalent to 

^+tV=0 (Bll) 

P 

'Since //and s are constant along stroamlines (ref. 10), this approximates the general 
small-perturbation flow to the extent that the lines y=conslant approximate streamlines. 


which is just the linearized Bernoulli equation in terms u' 
alone. 

The physical interpretation is this: If the assumed total per- 
turbation consists oj a plane shear flow (u”, 0) and a potential 
flow ill' , r'), then there is no pressure perturbation associated 
with the shear flow; the entire pressure perturbation arises from 
the potential flow and is related to u' by the ordinary linearized 
Bernoulli equation. In other ivords, the pressure is obtained by 
subtracting out the shear-flow velocity and applying the linear- 
ized Bernoulli equation to the remaining velocity. 

DENSITY FIELD 

The density distribution can be related to the velocity and 
entropy distributions as follows: fl'lie starting point is the 
dift’erential equation of state (B2) 

5 p bp 5s 

P pa'^ Cp 

Again assume that the flow is a combined potential flow 
(u',v') and shear flow (See eqs. (BlO) and after.) 

Then equation (Bll) applies for bpjp, and the density field 
is given by 

5p_ Uu' 5s 

P~ a^ Cp 

or 

(B12) 

ThuJt it is found that the density perturbation depends on the 
potential flow via the velocity perturbation u' and on the shear 
flow via the entropy perturbation 5s. 


« Since II and s are constant along streamlines (ref. 10), this approximates the general 
sniall-perlurbation flow to the extent that the lines ii=constant approximate streamlines. 


APPENDIX C 


VARIANTS OF THE SHOCK RELATIONS 


The ratio of the normal velocities before and after the 
shock has been defined as m: 


m = UJU 

Thus by reference 8, equation (114), 

7+1 y; 2 


(Cl) 


(C2) 


where Ua is the normal Mach number ahead of the shock. 
Corresponding!}', U is the normal Mach number behind the 
shock, and by reference 8, equation (112), 






77 2 7—1 

70 A s 


(C3) 


From equations (C2) and (C3) it can lie shown that 

U'- 2 


and 


clier 


(7+l)(m-l) 

r+-i 




— 1 


(C4) 

(C5) 




Thc equalit.v of transvei-se velocit.v components across an 
oblique-shock wave requires, in the present notation, that 

Va tan d=l~ tan ip 


Then, with the definition (Cl), 

m tan 0=tan a (C6) 

Equations (C2) and (C6) together allow ip to be determined 
in terms of d and Ua- 
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